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Abstract 

We studied the phononic heat transfer through an atomic dielectric wire 
with both infinite and finite lengths by using a model Hamiltonian approach. 
At low temperature under ballistic transport, the thermal conductance con- 
tributed by each phonon branch of a uniform and harmonic chain cannot 
exceed the well-known value which depends linearly on temperature but is 
material independent. We predict that this ballistic thermal conductance will 
exhibit stepwise behavior as a function of temperature. By performing nu- 
merical calculations on a more realistic system, where a small atomic chain is 
placed between two reservoirs, we also found resonance modes, which should 
also lead to the stepwise behavior in the thermal conductance. 
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I. INTRODUCTION 



Many studies have demonstrated that transport properties of condensed matter with 
nanometer size range and low dimensionahty can exhibit interesting features. For example, 
the electron transport through a nanoobject (i.e. a small molecule, nanowire or atomic chain) 
shows novel feature JiHl owing to two facts: (1) The electronic states confined in directions 
which are perpendicular to the one for the current flow are quantized. The level spacing^, 
Aej = ej+i — ej, of the quantized states of the nanoparticle are finite but strongly size and 
geometry dependent. (2) The mean distance, /, traveled by an electron between two termi- 
nals is smaller than the phase breaking length, L<^, (i.e. / < L^). Consequently, a stepwise 
behaviour in the electrical conductance is observed under ballistic transport through a point 
contact created by an STM tipil0 (or through a narrow constriction between two reservoirs 
of two-dimensional electron gas in high mobility GaAs-GaAlAs heterostructureslli), because 
the quantization and coherence of electronic states become pronounced. The ballistic elec- 
trical conductance G was foundHi to vary with the radius of contact, rc, (or with the width 
of the constriction, w) in quantized steps of approximately 2e^ /h, when or w is in the 
range of Fermi wave length Xp- 



In view of the interesting physics underlying the ballistic electron transport through a 
nanoobject between two electrodes, we conjectured! earlier that similar effects should occur 
also in the heat conduction via phonons. Similar to a constriction in electrical conduction, 
we expect that a nanoobject has a discrete vibrational frequency spectrum with finite spac- 
ings, AQi = Qi+i — Qi, determined by its atomic configuration and size. Furthermore, when 
a nanoobject is placed between two reservoirs, collective vibrational modes (phonons) which 
involve the coherent vibrations of the atoms in the nanoobject and reservoirs can be devel- 
oped, leading to the process of ballistic thermal conductivity. The thermal conductance, /C, 
depends on the available number of channels and hence it is expected to increase in discrete 
values each time a new ballistic channel is included in the heat transfer. Each vibrational 
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mode contributes to the heat transfer, but its contribution is weighted by Planck's distri- 
bution n{Q,T). Therefore, /C is determined by phonons with energies up to thermal energy 
{hQ ~ ksT). This may lead to smearing of the step structure of /C. This feature of quanti- 
zation of thermal conduction differs from that of electrical conduction, because the former 
complies with Bose-Einstein statistics. 

Heat transfer through a nanoobject is also related to the atomic scale dry friction and 
lubrication0"Ei. When two objects move relative to each other, they engage in contact 
through various atomic scale asperities or lubricant molecules. In these studies, the em- 
phasis is placed on the energy dissipation. The mechanical energy of the relative motion 
is dissipated by (a) the creation of electron-hole pairs, and (b) by the creation of non- 
equilibrium phonon distribution at the close proximity of the asperityElS. In the latter 
case, the main issue is the thermal conductance related to the heat transfer between the 
vibrational modes of the asperities (or lubricant molecule) and those of the objects, that 
leads to the energy dissipation^. In addition, the heat transfer is also an interesting issue 
in nanodevicesEl. Therefore, it is important to understand the thermal conduction through 
a nanoobject. 

The objective of this work is to investigate the effects of the quantized nature of vibra- 
tional modes of a dielectric wire on phononic heat transfer between two reservoirs. To this 
end, the ballistic heat transfer through a uniform dielectric wire is investigated first. We 
also studied the thermal conduction through a finite atomic wire. As an example, the heat 
transfer through a three-atom chain placed between two reservoirs with large cross section 
are examined by performing numerical calculations on a model system^. We found that 
the thermal conductance of a uniform atomic wire at low temperature is independent of 
material and increases linearly with temperature. For finite atomic wire both longitudinal 
acoustic and optical modes contribute to the thermal conductance; the latter modes may 
lead to stepwise behavior. 

3 



II. MODEL HAMILTONIAN APPROACH 



The system of interest is schematically described in the upper panel of Fig. 1. Heat flows 
from the left (L) to the right (R) through a nanoobject (C). Here L, R can be sections of 
nanoobject itself or heat reservoirs. The temperature of L and R, Tl and Tr, respectively 
are kept constant, so that the steady state of the transport is considered. For this system 
we investigate the thermal conduction by using a model Hamiltonian approach. To find the 
modes which are responsible for the ballistic transport of heat, we specify first a zero order 
vibrational Hamiltonian which characterizes the possible modes in the three regions, (L + 



Here qx is the momentum quantum numbers for vibrational mode ^iqx belonging to the 
branch i in one of the three regions (X = L,R,C), and Ojg^ i^tqx) corresponding 
annihilation (creation) operator. The heat transfer via phonons is realized when L,C, and 
R are coupled and AT = — Tr > 0. The coupling can be described by an interaction 
Hamiltonian, 



ijt 

in terms of the harmonic (T) and anharmonic (U) coupling coefficients. This Hamiltonian 
has been used to predict localized and extended vibrational modes in molecular crystals^!. 
The present application is similar to the transfer Hamiltonian approachlll. As we empha- 
sized earlier, we consider only the thermal conduction mediated by phonons. The electronic 
states, as well as the electron-phonon interactions are not included in Eq.(l) and Eq.(2). 



C + R), 



Ho — ^ ^^iqL^tqL^iqL 5Z ^^iQR^tqR^iqR ~^ ^ ^^iqc^tqc^^lC 



(1) 





(2) 



The modes responsible for the ballistic heat transport are provided mainly by the har- 
monic couplings. Before discussing the effects of the harmonic coupling, we comment on 
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the anharmonic terms. The anharmonic couphngs describe the phonon-phonon interactions 
which are the mechanism responsible for several dissipative processes, such as energy loss 
to local modes of the nanowire and thermal resistance due to umklapp process. Depending 
on the strength of the anharmonic coupling, the continuing energy transfer to the local- 
ized modes may result in energy discharge to the heat transporting states, or desorption 
or migration of an atom. The former event, i.e. the excitation of a local phonon, followed 
by deexcitation can occur periodically and can be manifested in the time variation of the 
conductance, /C(t). The detail of the anharmonic effects will not be discussed further in this 
paper. In what follows we apply the above approach to calculate the thermal conductance 
of a uniform dielectric wire (Sec. Ill) and also of a finite atomic chain (Sec. IV). 



III. UNIFORM DIELECTRIC WIRE 

We now consider three sections on a uniform wire. Owing to the uniformity of the sections 
(L, C, R), the phonon modes are the normal modes of the system, hence q = Ql = Qc = Qr, 
and also Tij^gg = Ta^qq. Furthermore, %i,qq = hQiq. If the effect of mode decay is considered, 
the rate of the heat transfer from L to i? in this ID system is given by Jq^ = J2iq ^^iq^iq 
by using the decay rate expressionEl for a given i and q, 



yLCyCL 

■pL _ pC qq qq (o\ 



with rg = i2n/h)Vq^/Vq^^Dj'{Q,q). By symmetry V^f^ = V^f and 

Ka*^ = {r^L -l,nc + 1, nnlHintlriL, nc, ur) 



Vgg = {riL - 1, nc, ur + l\HM\nL -l,nc + 1, hr). (4) 

where nx {X = L,C,R) is the Planck's distribution nx{^iq,Tx) = [exp{hQiq/ kBTx) — 1]^^ 
given in terms of the frequency of the mode and temperature of the section, Tx- Hence 



Vq^q'^ = %i,qq\JnL{nc + 1) aud V^q^ = Tii^qq^J{nc + l)inR + l) are expressed in terms of 
Planck's distribution functions. The net rate of heat transfer through the chain is Jq = 
Jq^ — Jq^ ■ Since and n/j <^ 1 at low temperature, Jq is expressed in the following form 



Jq-ttT. [nL-nR]ndn (5) 

ZTT ^ JO 

Here index i of the summation sign specifies the branches for ul and n/j. Note that 
Eq.(5) yields exactly the same expression of current one obtains from the phenomenological 
approachiiEll by assuming the transmission coefficients, ti{fl) = 1. Therefore, the present 
approach justifies the Landauer type expression of heat current by deriving it from the model 
Hamiltonian specified by Eq.(l). Now, we consider first an atomic chain that has only three 
branches and evaluate Eq.(5) for AT 0. We take Tl and Tr small so that the upper limit 
of the integral in Eq.(5) is justified. Defining T = (Tl — Tr)/2, the following expression for 
the ballistic thermal conductance, /C = Jq/AT, is obtained, 

(6) 

According to this interesting result the ballistic conductance of each branch i of the uniform 
and harmonic atomic chain is limited by the value JCq = Ti^k\T /"ih. It is independent of any 
material property, and is linearly dependent on T. The total thermal conductance becomes 
K, = MlCoi where M is the total number of branches. For an ideal ID atomic chain A/" = 3, 
if the transverse vibrations are allowed. The result agrees with the well-known expression 
derived earlier for dielectric bridges&^. For a homogeneous wire that has a larger cross sec- 
tion, more subbands can now be present due to the transverse confinement of the vibrational 
motions in the direction perpendicular to the propagation. In this case, the ballistic thermal 
conductance of a branch with VLi^min 7^ is smaller than /Co due to the finite frequency of 
a subband through the occupancy of Planck's distribution. For example, with two subband 
frequencies, VLmin = lO^^Hz and Qmin = lO^^Hz at q = 0, the thermal conductance of each 
channel at T = 0.25K is calculated from Eq.(5) to be 0.94/Co and 0.25/Co, respectively. At 
T = IK, these contributions increase to 0.97/Co and 0.77/Co, respectively. Consequently, one 
expects that each subband in the wire having Qi^min 7^ 0, that formed by increasing the cross 
section, leads to a jump in /C as in the electrical conduction. However, the amplitude of steps 
decreases exponentially with increasing Qi^min at given low temperature T. Furthermore, an 
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important difference between the electrical and thermal conductance is that all the modes 
of a given branch i contribute to /Co, whereas the quantum ballistic electrical conduction in 
the ID channelsi'B is governed by the Fermi-Dirac distribution. 

Recently, Greiner et al.il calculated the thermal conductance of free electrons for ID bal- 
listic and diffusive transport in the linear response regime by using the correlation function 
formalism. They derived a universal expression for the electronic thermal conductance of a 
degenerate electron gas, Kf, that is the same as /Co found for the phononic thermal conduc- 
tance in Eq.(6). Interestingly, electrons and phonons have the same universal value for their 
thermal conductance at low temperature, in spite of the different statistics they obey. We 
believe that this situation originates from the fact that both systems follow the same distri- 
bution, e"^/'^^^ for e/ksT >> 1. Here, e = h^li for phonons, and e = ti — fi for electrons, 
and /i being electron energy and chemical potential. 

IV. FINITE WIRE AND NUMERICAL MODELING 

The quantum heat transfer through a uniform and harmonic atomic wire is an idealized 
situation. A more realistic system consists of a finite atomic wire, C coupled to reservoirs 
L and R. Here the finite atomic wire has discrete vibrational frequency spectrum, while 
those of reservoirs are continuous. In the lower panel of Fig. 1 the corresponding densities 
of states, D-^iyi), are schematically described. It is also noted that cross sections change 
abruptly where the atomic wire is connected to the reservoirs. As a result, the transmission 
coefficient, tj{VL) < 1 because of reflections. On the other hand, the application of the con- 
tinuum theory to determine the transmission coefficients as in the case of dielectric wiresH 
having cross section (50nm x 50nm) is not justified here owing to the discrete character and 
small size of the finite atomic wire. Moreover, the modes of wire differ from those of the 
reservoirs; C itself is not uniform in the sense that the atoms near the boundary experience 
different force constants. We now treat such a general non-uniform dielectric wire. 
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In order to analyze the character of the modes of a finite wire and to present an exam- 
ple we carried out numerical calculations on a prototype model that consists of a chain of 
3 atoms placed between two reservoirs. Each reservoir (L and R) is modeled by 7 (001) 
atomic planes of the fee structure. The whole system (L + C + R) comprises 451 atoms. The 
interatomic distances between an atom in the chain and nearest slab-atoms and between 
its adjacent chain atoms are taken, 2.7, 2.6A, respectively. We assume that the electrical 
conduction through these atoms is zero, and hence the energy transfer through electrons is 
not allowed. Since we are concerned with the vibrational modes of the system, which are 
responsible for the energy transfer, we express the interatomic interactions by a prototype 
empirical, Morse-type two-body potential. Then, we set up the dynamical matrix of the 
system and calculate the normal modes and their frequencies flig. 

Our analysis distinguished three different types of modes: i) Modes localized strongly in 
the chain, ii) Modes confined to the reservoir (s). iii) Current-transporting resonant modes. 
Note that only confined modes would appear if the coupling between the chain and the 
reservoirs were weak. The modes in the first two categories ( i, ii) cannot be involved in the 
heat transfer, but they affect the thermal conduction by increasing the thermal resistance 
through anharmonic couplings. The transverse modes of the finite chain are well localized 
owing to the different coupling constants and different number of neighbors as compared 
to the atoms in the reservoirs. On the other hand, the longitudinal modes are coupled to 
those of the reservoirs and form resonant modes. Figure 2 illustrates three resonant modes 
of the atoms. Their relative motions are indicated by the arrows, (a) and (b) can be consid- 
ered as longitudinal optical branches and (c) the acoustic branch. Therefore, the resonant 
modes can be also contributed by the optical branches. Figure 2 also shows that the local 
densities of states within the chain, pjg^(Q), associated with these modes, iqc, which are 
calculated from the spectrum of the dynamical matrix. Wc now consider the ballistic heat 
transfer through resonances. Starting from the rate expression given in Sec. Ill, the ballistic 
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thermal conductance for a single resonant mode iq, is given by, 

Kiq ~ kBViqx'^e''(e'' - 1)'^ (7) 

with X — hO^iq/kBT. For large x, Kiq ~ Viqksx'^e'^ . On the other hand Kiq ~ viqkB if x is 
small. In view of the analysis in Ref [25], the thermal conductance of a finite atomic chain 
due to all resonance modes can be obtained from the sum over these modes, iqc-, 

-2 

dVt (8) 

Here piq^ is the local density of states for the resonant modes derived from the mode iqc 
and Viqc is average velocity associated with the same mode. This density of states modifies 
the ideal ballistic transmission as tiq{Q) does to Jq in Ref [23]. If p^g^, in Eq.(8) is replaced 
by Dfiyt) corresponding to a branch of uniform atomic chain, the thermal conductance AT 
in Eq.(8) yields the universal value /Cq. By using a value of Viq^ ~ 393m/sec deduced from 
the slab calculations, we calculate /C(T) due to resonance in Fig. 2(c) and illustrate it in 
inset in the same panel. For T << 50 K, /C increases with T, then /C saturates as T — > 50 
K. This latter fact is the consequence of Kiq at small x discussed above. 
We shall emphasize an interesting feature of /C that can be revealed from Eq.(8). Each atom 
incorporated in the finite chain between two reservoirs adds 3 (optical) modes. So as the 
chain becomes longer by including new atoms, the number of longitudinal optical modes 
which may give rise to resonant modes increases. Accordingly, each resonant mode labeled 
by iqc that can lead to a current transporting channel, shall give rise to a step increase in 
K, by increasing the temperature. Apparently, the step behavior can occur without having 
subbands in the transverse direction. This situation is relevant for the energy transport 
through molecules attached (i.e chemisorbed or physisorbed) to the surfaces. The ther- 
mal conductance associated with the phononic heat transfer for a molecule having discrete 
vibrational frequency spectrum may exhibit at least non-monotonic variation with temper- 
ature. Whether the step behavior of thermal conductance can be resolved depends on the 
frequency and level spacing of the resonance modes. According to a criterion put forward 
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by BlencowetJ, Qi+i > lAQi in order to resolve the quantum features of /C(T). As for the 
question whether /C of a molecule or an atomic chain between two reservoirs increases by 
incorporating new atoms, the answer is not straightforward. It depends on the new modes 
yielding new resonance mode, and also piq^{Q) and Viq^ associated these modes. Neverthe- 
less, this is an interesting issue that deserves further study. The contacts to the reservoirs 
by themselves can set a resistance. It should also be noted that increasing the length of 
the chain by including more atoms can lead to the possibility of introducing imperfections 
in the structure. The phonon-phonon interaction due to anharmonic coupling (that is not 
considered here) can also contribute to the resistance of the wire. Hence the probability of 
observing such coherent conductance decreases. As a result, the resistance to heat current 
increases. 



V. CONCLUSIONS 

In this work we analyzed important features of the heat conduction via phonons through 
a dielectric wire. In particular, we investigated the ballistic heat transfer through a uniform 
and harmonic wire and also through a finite atomic chain. By using a model Hamiltonian 
approach we showed that the conductance of a uniform atomic chain at low temperature is 
consistent with the well-known universal value, and varies linearly with temperature. As the 
cross section of the wire becomes larger to include several subbands due to the confinement of 
the transverse motion of ions, each subband contributes to the ballistic thermal conductance 
by an amount that depends on T and the threshold frequency of the subband, fii,mm- The 
higher the frequency of the mode, the less its contribution to /C. The resulting step behavior 
shall be observable, if they are not smeared out. We analyzed also the nature of modes for 
a finite atomic chain placed between two reservoirs. We found resonance modes originating 
from longitudinal optical modes, which can contribute to the ballistic thermal conductance. 
They differ from the subbands due to lateral confinement but should also exhibit stepwise 
behavior in /C. 
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FIGURES 

FIG. 1. A schematic description of the quantum heat transfer investigated in the text. A 

nanoobject coupled to two reservoirs and corresponding density of frequencies. Resonance states 
are described by dotted hues. 

FIG. 2. The local density states corresponding to resonance modes derived from the three 
longitudinal modes iqc^s of a chain of three atoms between two slabs representing the reservoirs 
are illustrated in panels (a)-(c). The corresponding modes are schematically described at right. In 
(c) the variation of /C in units of 10~^^ Jm/K sec calculated from the acoustic resonance modes is 
illustrated in the same panel. 
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